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Abstract 

We study the contributions So and Si, proportional to a and a , to the fermion self- 
energy in Wilson's formulation of lattice QCD with UV-filtering in the fermion action. 
We derive results for m cr jt and the renormalization factors Zs, Zp, Zy, Za to 1-loop order 
£NJ ■ in perturbation theory for several filtering recipes (APE, HYP, EXP, HEX), both with 

and without a clover term. The perturbative series is much better behaved with filtering, 
in particular tadpole resummation proves irrelevant. Our non-perturbative data for m cr it 
and ZA/{Z m Zp) show that the combination of filtering and clover improvement efficiently 



O 



reduces the amount of chiral symmetry breaking - we find residual masses am rcs = 0(1CP 2 ) 

> 

1 Introduction 

The Wilson formulation of lattice QCD breaks the chiral symmetry among the light flavors pD,[2]. 
Accordingly, Wilson fermions undergo an additive (and multiplicative) mass renormalization. 
While this is not a problem in principle - the explicit breaking disappears if the lattice spacing a is 
sent to zero [3] it entails a number of complications in numerical work based on this formulation. 
There are several strategies how the additive mass renormalization might be reduced. A popular 
choice, to augment the action by a clover term, has the merit of reducing cut-off effects from 
0(a) to 0(ag 2 , a 2 ) [U El E]. Another possibility, referred to as UV-filtering, is to replace all 
covariant derivatives in the fermion action by smeared descendents, as proposed in a staggered 
context El El and later applied to Wilson/clover fermions PHI HH El EH E] . We find that 
filtering indeed ameliorates important technical properties of the Wilson operator, as does the 
clover term without filtering. The real improvement, however, comes from combining the two. 
With standard conventions the (r = l) Wilson operator takes the form 

Dw(x, y) = - {(7^ - ^)^(>)<W, a - (7 M + I)Ul(x - + ir-8x, v (1) 

where I is the identity in spinor space. The Sheikholeslami-Wohlert "clover" operator follows 
by adding a hermitean contribution proportional to the gauge field strength [5] 

£>sw(z, y) = D w (x, y) - 8x, y (2) 

with = 7=jh7t,7*,] and F^ u the hermitean "clover-leaf" operator. In order to cancel the 0(a) 
contributions, the coefficient Csw needs to be properly tuned. In perturbation theory one finds 
csw = 1 at the tree-level and a correction proportional to the n-th power of at the n-loop 
level. It is well known that for the standard "thin link" operator perturbation theory shows 
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thin link 1 APE 2 APE 3 APE 1 HYP 


c sw = 
csw — J- 
csw = 2 


51.43471 13.55850 7.18428 4.81189 6.97653 
o±.yo044 4.yuoio 1.004O0 u.Muyo i.yoooi 
1.10790 -7.11767 -5.48627 -4.23049 -4.41059 



Table 1: Additive mass shift 5* for "thin link" Wilson or clover fermions and after APE or HYP 
filtering with standard parameters. The uncertainty is of order one in the last digit quoted. 



rather bad convergence properties. Therefore, the ALPHA collaboration has started a non- 
perturbative improvement program [15]. Another approach is to resum the tadpole contributions 
|16j . since they are quite sizable. For filtered Wilson/clover quarks this might be different - we 
elaborate on "fat link" perturbation theory (TTJ [171 OS EH], and we compare these predictions 
to (non-perturbative) data. It turns out that filtered perturbation theory shows a much better 
convergence behavior, but still, it does not describe the data very accurately. The agreement 
is (at accessible couplings) much better than in the unfiltered theory, but it is far from being 
completely satisfactory. We find that the additive mass shift is two orders of magnitude smaller 
than without filtering, and this is extremely useful in phenomenological studies. 

The following two sections contain our perturbative results for UV-filtered Wilson/clover 
fermions. Sect. 2 focuses on the additive mass shift with 1,2,3 steps of APE, HYP, EXP, HEX 
filtering and arbitrary improvement coefficient csw- Sect. 3 contains our 1-loop results for the 
renormalization factors Z$, Zp, Zy, Za with these filterings, a reminder how improved currents 
are constructed, and a comment on tadpole resummation. Sect. 4 presents our non-perturbative 
data for the additive mass shift and some renormalization factors, both with csw = and csw — 1- 
Sect. 5 contains our summary. Details of "fat-link" perturbation theory, of an explicit mass shift 
calculation and of the parameter dependence have been arranged in three appendices. 



2 Additive mass shift with UV filtering in 1-loop PT 

In this paper we consider four types of filtering: APE, HYP, EXP, HEX. The fist two are 
well known [19[ H5J, the third one has been named "stout" in [20], and the fourth one is a 
straightforward application of the hypercubic nesting trick on the latter (see App. A for details). 
While on a technical level the smearing produces a smoothed gauge background, it is in fact 
a different choice of the discretization of the covariant derivative in the Dirac operator and 
therefore leads to an irrelevant change of the fermionic action (provided the filtering recipe is 
unchanged when taking the continuum limit). 

In our analytical and numerical investigations we use the "standard" parameters 

^APE n R ^EXP n i / Q \ 

a 8 td = - 6 > «std = - 1 ( 3 ) 

for APE and EXP smearing, and similarly the "standard" parameters 

a s H J p = (0.75, 0.6, 0.3) , a™ x = (0.125, 0.15, 0.15) (4) 

for HYP and HEX smearing. The two values in ([3]) are related by giving an identical 1-loop 
prediction for all quantities of interest (e.g. — am crit ), and the same statement holds for the 
hypercubically nested recipes (T4|), see App. A for details. Accordingly, all perturbative tables 
with label "APE" will apply to EXP, too, and ditto for a label "HYP" and the HEX recipe. 
The additive mass shift is given by the self-energy S via [note that am CT \t < with ([TBI ] 

am crit = S = S~ 2 C F S + O(g^) (5) 

1671^ 
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where S is the quantity that is usually tabulated and Cj? = 4/3 for SU(3) gauge group. Gener- 
alizing a standard calculation [21] to "fat-link" perturbation theory (see App. A for a summary) 
one may work out 1-loop predictions for S [HI [TS]. We have done this for arbitrary cgw From 
inspecting Tab. [TJ one notices that csw = l alone reduces the additive mass shift by a factor 1.6. 
Filtering alone achieves a factor 3.8 or 7.4 with a single APE or HYP step, respectively. However, 
the combination reduces it by a factor 10.5 or 26.0, and hence proves much more efficient than 
any one of the ingredients alone. The tuned csw that would achieve zero mass shift is slightly 
above 2 in the thin-link case, and slightly above 1 in all cases with filtering. This is the first 
indication that filtered cgw = 1 clover fermions break the chiral symmetry in a much milder way 
than filtered Wilson or unfiltered clover fermions. An important question is, of course, to which 
extent this is realized non-perturbatively, and we shall address this issue in due course. 



3 Renormalization factors with UV filtering in 1-loop PT 



3.1 Generic setup 

In general, the matrix elements of some operator 0™ nt (/i) in the continuum MS scheme and its 

la1 
k 



lattice counterparts 0^. att (a) are related by 



(•|or(/i)l-) = E^M(-lol att (a)|.) 



(6) 



Zjk(afi) = 5 



9l 



(A 



latt 



A cont \ 

with Cp =4/3 for SU(3) gauge group. Typically (e.g. for 4-fermion operators and a non-chiral 
action), k runs over other chiralities than j. For 2-fermion operators, this mixing shows up at 
higher orders in an expansion in the lattice spacing a, and packing it into the construction of 
improved currents, one is left with the diagonal term in 0. With our convention (which agrees 
with [21 J , but not with [14J) a value zx>0 signals Z X <1- Specifically (with X — S, P, V, A), 



9l 



(7) 



Zp(afi) 



1 - 



9l 



4vr 2 



z s 



- log(a V 



1 - 



9o \ Z P i / 2 2\ 



12tt 2 
12tt 2 



zv 



za 



(9) 



for the (pseudo-) scalar densities and the (axial-) vector currents, with corrections of order 0(qq ) 
throughout. 



3.2 Results for Zs, Z P , Z v , Z^ for Wilson and clover fermions 

The same approach of combining FORM-based [22] standard perturbative procedures [21] with 
"fat-link" perturbation theory that has been used in the previous section for the additive mass 
shift, allows one to work out the renormalization factors Zs, Zp, Zy, Za for arbitrary csw- 

Our results for zx with X = S, P, V, A in the unimproved case csw — are summarized in 
Tab.[2J An important check is that (zp — zs)/2 and zy — za should coincide [23J. The pertinent 
entries indicate that the integration routine yields at least 6 significant digits. 

Our results for zx with X = S, P,V, A in the improved case csw = 1 are summarized in Tab.[3j 
Again we check the quality of the agreement between {zp — Zs)/2 and zy — za- Moreover, since 
these figures indicate the amount of chiral symmetry breaking [23], it is instructive to compare 
the bottom lines of Tab. [2] to those of Tab.O Improvement alone reduces zy— za by a factor 3.1. 
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c sw = 


trim link 


1 APE 


2 APE 


3 APE 


1 TT\7""n 

1 HYP 


zs 


12.95241 


1.12593 


-1.53149 


-2.87223 


-1.78317 


zp 


22.59544 


5.28288 


1.07019 


-0.98025 


0.51727 


zv 


zU.ol (ok) 


o.oyoiu 


o.ozzol 


z.olool 


o.ooUf 


za 


15.79628 


4.31963 


2.32197 


1.56782 


2.23054 


\zp-z s )l I 


/i soi c;o 


O H78 /I S 


1 Qnns/i 
l.oUUo4 


U.y4oyy 


1 1 KflT) 

1. loUzz 


zv-za 


4.82152 


2.07847 


1.30084 


0.94599 


1.15022 



Table 2: Coefficient ,2x in formula (j7]) for the renormalization factor Zx with X 
c sw = Wilson fermions with APE or HYP filtering with standard parameters. 



S, P, V, A for 



csw = l 


thin link 


1 APE 


2 APE 


3 APE 


1HYP 


1 HYP [14J 


zs 


19.30995 


4.11106 


0.40606 


-1.43930 


-0.03678 


0.12 


zp 


22.38259 


4.80364 


0.65185 


-1.33218 


0.12845 


-0.04 


zv 


15.32907 


3.31243 


1.43934 


0.82550 


1.38517 


1.38 


za 


13.79274 


2.96614 


1.31645 


0.77195 


1.30255 


1.30 


(z P -z s )/2 


1.53632 


0.34629 


0.12290 


0.05356 


0.08262 


-0.08 


z v ~z A 


1.53633 


0.34629 


0.12289 


0.05355 


0.08262 


0.08 



Table 3: Like Tab.[2l but for csw = l clover fermions. The last column has been adapted to our 
sign convention [cf. (J7|)] and suggests that there is a mislabeling in Tab. Ill of Ref. EE 



csw = 2 


thin link 


1APE 


2 APE 


3 APE 


1HYP 


zs 


22.90672 


4.35133 


0.06571 


-1.91937 


-0.43671 


zp 


26.24177 


6.10928 


1.39146 


-0.81914 


0.80287 


zv 


8.95400 


-0.33664 


-1.07948 


-1.08366 


-0.89073 


za 


7.28648 


-1.21561 


-1.74236 


-1.63378 


-1.51052 


(z P -z s )/2 


1.66753 


0.87898 


0.66288 


0.55012 


0.61979 


z v -z A 


1.66752 


0.87897 


0.66288 


0.55012 


0.61979 



Table 4: Like Tab. [2], but for csw = 2. This nails down the full polynomial dependence on csw- 



One step of APE or HYP filtering diminishes it by a factor 2.3 or 4.2, respectively. However, the 
combination of these recipes achieves a factor 13.9 or 58.4, and hence proves much more efficient 
that any of the ingredients alone. This is in line with the lesson learned from Tab.[U 

Our results for zx in the case csw = 2 are shown in Tab. HI Obviously, "too much" im- 
provement deteriorates the chiral properties of the action. At 1-loop order all zx depend on 
Csw through a quadratic polynomial, hence Tabs.[2]-H] give them for arbitrary values of the 
Sheikholeslami-Wohlert parameter. For instance, for 1 HYP (or 1 HEX) step they imply 

Zs = -1.78317 + 2.81955c sw - 1.07316c| w 

Zp = +0.51727 - 0.92044c sw + 0.53162ci w 

z v = +3.38076 - 1.85544c sw - 0.14015c| w 

z A = +2.23054 + 0.01455c sw - 0.94254c| w (10) 

and from the pertinent curves (see Fig.[T]) one learns two lessons. First, the point where the 
1 HYP action is most chiral (i.e. where zp — zs and zv — za are minimal) is near csw — 1-1653. 
Second, near c-sw = l-5 the four coefficients zs,p t v,A are simultaneously small. By contrast, with 
less filtering (e.g. 1 APE) the point of minimal chiral symmetry breaking is further away from 
1, and the four renormalization factors cannot be simultaneously close to 1. 
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0.5 1 1.5 2 0.5 1 1.5 2 

c sw c sw 

Figure 1: Finite pieces z^py^ of the Zx for 1 APE and 1HYP fermions as a function of csw- 

Any strategy in which csw deviates, for large j3, from 1 by a polynomial in g$ with vanishing 
constant part yields a theory with 0(agl) cut-off effects. Here we restrict ourselves to csw = 1- 
Getting higher terms in the polynomial right reduces discretization effects to 0(ag§) or better, 
and non-perturbative improvement would realize 0(a 2 ). 



3.3 Construction of improved currents and densities 

At tree-level Zs,py,A — 1> an d the improvement coefficients are csw — 1> ^s,p,v,a — 1> b m = —1/2 
and cy jV 4 = 0. Accordingly, in a tree- level 0(a) improved theory the currents read 

(S imp ) a = (l + am q )S a 

(P imp ) a = (l + am q )P a 

(V imp )« = (l + am q )V* 

(Amp); = (l + am s )A; (11) 

which is free of mixing effects, but it is well known that (at least in the unfiltered case) this is 
not sufficient to be in the Symanzik 0(a 2 ) scaling regime for accessible couplings. Throughout, 
we use the flavor decomposition X = X a ^- with A a one of the Gell-Mann matrices (a=1..8). 

At the 1-loop level and with Nf = 0, N c = 3, renormalization factors in the unfiltered theory^ 
take the fornix = 1-0.163042^, Z P = 1-0.188986^, Z v = 1-0.129430^, Z A = 1-0.116458^, as 
follows from the first column of Tab.H Similarly c sw = 1 + 0.2659c^ |2l] , and b s = 1+0.1925^, 
b P = 1 + 0.1531^, b v = 1 + 0.1532^, b A = 1 + 0.1522^, b m = -1/2-0.09623^, c v = -0.01633^, 
ca = — 0.00757(7q, see [151 [25], [26] [27] for details. The main message is that most of the 1-loop 
corrections are large, since 1. With these expressions at hand, improved currents follow via 



('S'imp 


r 


= Z s ~S a , 


s a = 


(l + b s am q ) 


S a 


(-Pimp 


r 


= z P P a , 


pa 


(1 + bp am q ) 


pa 


(Mmp, 


° 

A" 


= ZyV; , 


v a = 


(1 + by am q ] 


[v; + ac v B u T; v 


(^imp, 


Q 
/' 


= z a ai , 




(1 + b A am q ) 


[Al + ac A d,P a ] 



(12) 

where <9 M = |(<9 M + <9*) denotes the forward-backward symmetric derivative. Clearly, this is a 
complicated mixing pattern involving even the tensor current. Still, with perturbative coefficients 



1 Throughout, we use csw to the previous order in quantities which depend on it; these Zx are for csw = 1- 
2 With Nf > they depend on <jf§ = + b g am w ) with b g = 0.012000(2)^/ and m w given in (fT3 j) [IS ] . 
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it remains (in the unfiltered theory) a challenge to reach those couplings where the Symanzik 
scaling with 0(a 2 ) cut-off effects sets in. This is why (in the thin-link theory) a non-perturbative 
determination of the renormalization constants and improvement coefficients is preferred |15j . 

Our hope is that with filtering perturbative improvement at the 1-loop level is a viable 
strategy. An important check is how well the renormalized VWI quark mass and the renormalized 
AWI quark mass coincide. The (bare) Wilson or clover quark mass is defined as 

m w = m - m crit where am = \ (- —) , am crit = \(— —) (13) 



tree 



„ 1 I j UU „ 1 

Z K Ktree ^ ^crit K tr 

with K tree = l/8, and the (renormalized) VWI quark mass then follows through 

m VWI (/x) = Z m (a/i)(l + b m am w )m w . (14) 

The (bare) PC AC quark mass is defined through (for and P built from degenerate quarks) 

pcac _ l (d,[Al{x) + ac A d,P a }O a m 
m "2 <P*(x)O*(0)> 1 ' 

and the (renormalized) AWI quark mass then follows through 

m AWi/ s = Z a 1 + b A am w PCAC 

In (fT4l [T6|) the details of the conversion from the specific cut-off scheme on the r.h.s. to the 
standard MS-scheme on the l.h.s. are built into the renormalization factors. If we had csw and 
the bs, bp, by, b A , b m , cy, c A at 1-loop level, plus the Zg, Zp, Zy, Z A at 2-loop level, a theory with 
0{ago) cut-off effects could be realized. At the time, we lack the knowledge of any improvement 
coefficient at the 1-loop level (with filtering). Accordingly, the following section is devoted to 
a preliminary test with tree-level improvement coefficients and 1-loop renormalization factors. 
Still, since the perturbative series converges so well, our hope is that this test does not fail 
completely - otherwise higher order corrections could barely save the case. 



3.4 Irrelevance of tadpole resummation 

One of the attractive features of filtered Dirac operators is that 1-loop renormalization factors 
and improvement coefficients are much closer to their tree-level values, suggesting a better con- 
vergence pattern. Obviously, a first guess says this is mostly due to the tadpole contribution 
being much smaller than in the unfiltered theory. 

In Feynman gauge the "thin-link" tadpole diagram with the value 12.233050gQCV/(167r 2 ), 
which is responsible for many of the large corrections in unfiltered perturbation theory [16], gets 
reduced as detailed in Tab. [5] for a broad range of a APE and n iter parameters. Note that these 
numbers hold for arbitrary csw, since the dependence on the Sheikholeslami-Wohlert parameter 
comes through quark-gluon vertices with an odd number of gluons. 

In Landau gauge the effect is even more pronounced, as shown in Tab.[6j Here, the "thin 
link" value is 9.174788(7qCf/(167t 2 ), and a smearing parameter a AFE <a^^ = 0.75 seems to be 
beneficial (cf. App. A for details on ct^ax)- I n this gauge the sunset diagram is rather small, 
regardless of the filtering level. We checked that, for the extreme choice (ct APE , 7iit er ) = (0.45, 10), 
we reproduce the result 0.2597053^C F /(16vr 2 ) of p]. 

From this observation it is plausible that tadpole improvement is not necessary - i.e. has 
barely an effect - in fat-link perturbation theory. This leaves us optimistic that the perturbative 
series might converge much better for filtered actions. The real issue is, of course, whether such 
perturbative predictions will agree with non-perturbative data. 
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0.12 


0.24 


0.36 


0.48 


0.6 


0.72 


0.84 


0.96 


1 


10.05384 


8.25363 


6.83240 


5.79017 


5.12693 


4.84269 


4.93744 


5.41118 


2 


c\ v r\c\ nr 

8.50285 


6.32137 


5.11011 


4.45066 


4.08470 


3.91401 


a r\ r\r\ a c\ 

4.00042 


4.56587 


3 


7.37921 


5.28658 


4.37748 


3.94939 


3.72277 


3.60854 


3.69503 


4.45440 


A 

4 


c rri A7 
D.OOlU 1 


A P.Q.AT7 


a nriQ/in 
4.UUo4U 


o 7n/i a i 
6. / U44 / 


o e.A 7 A A 
O.04 / 44 


o.4o/yo 


o.oo / zy 




5 


5.93054 


4.30964 


3.78626 


3.56346 


3.44535 


3.37845 


3.48725 


5.34886 



Table 5: Tadpole diagram in Feynman gauge [value to be multiplied with g^C f / (16n 2 )] in 1-loop 
"fat-link" perturbation theory. The corresponding "thin-link" value is 12.233050. 





0.12 


0.24 


0.36 


0.48 


0.6 


0.72 


0.84 


0.96 


1 


6.99558 


5.19536 


3.77414 


2.73191 


2.06867 


1.78443 


1.87918 


2.35292 


2 


5.44459 


3.26311 


2.05185 


1.39240 


1.02644 


0.85574 


0.94215 


1.50761 


3 


4.32095 


2.22832 


1.31922 


0.89113 


0.66450 


0.55028 


0.63677 


1.39614 


4 


3.49281 


1.62650 


0.94513 


0.64620 


0.48918 


0.40469 


0.49903 


1.64011 


5 


2.87228 


1.25138 


0.72799 


0.50519 


0.38709 


0.32019 


0.42898 


2.29060 



Table 6: Tadpole diagram in Landau gauge [value to be multiplied with g 2 Cp/ (16n 2 )} in 1-loop 
"fat-link" perturbation theory. The corresponding "thin-link" value is 9.174788. 



p 


5.846 


6.000 


6.136 


6.260 


6.373 


L/a 


12 


16 


20 


24 


28 


L/r 
a' 1 [ GeV] 


2.979 


2.981 


2.983 


2.981 


2.979 


1.590 


2.118 


2.646 


3.177 


3.709 


^-conf 


64 


32 


16 


8 


4 



Table 7: Matched (j3, L/a) combinations to achieve L/r = 2.98 as accurately as possible, based 
on the interpolation formula of [28J. n con { is the number of configurations per filtering and mass. 

4 Non-perturbative tests 

Here, we investigate how well a perturbative improvement program with 1-loop renormalization 
factors and tree- level improvement coefficients works with filtered Wilson/clover fermions. Since 
no phenomenological insight is attempted, we work in the quenched theory. We wish to cover 
a regime of couplings from (3 ~ 5.8 to /? ~ 6.4 with the Wilson (plaquette) action and we work 
in a fixed physical volume as defined through the Sommer radius r [28]. The corresponding 
parameters (realizing L/r = 2.98, and thus L~1.49fm if r = 0.5fm) are given in Tab. [71 

Technically, we produce a smeared copy of the actual gauge field, and evaluate the fermion 
action on that smoothed background. This differs from the approach taken in [T3], since our 
entire D w in (J2J) is constructed from smoothed links. See App. B for details. 

4.1 Data for m crit , Z A with APE/HYP/EXP/HEX filtering 

For clover fermions one has, up to 0(ag 2 , a 2 ) terms, the vector and axial- vector Ward identities 

Z v (.\d^{x)\.) = ZmM / sM (mf -<)(.|^(x+4) + 2S«(x) + S a (x-4)\.) (17) 
Z A (.\d^(x)\.) = ZmM / pM (m^+<)(.|P"(x+4) + 2P a (x) + P Q (x-4)|.) (18) 
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0.5 n 

APE1 i — I — 
APE3 h- * 
HYP1 i -»-- 
HYP3 : □ 
EXP1 
0-4 "EXP3 

HEX1 >---•-- 
HEX3 



0.5 n r- 

APE1 i — I — 
APE3 h- * 

HYP1 : ■-»-- 
HYP3 : □ 
EXP1 
0-4 -EXP3 
HEX1 : 
HEX3 



-0.25 -0.2 -0.15 -0.1 -0.05 0.05 0.1 0.15 




0.5 | 1 

APE1 ^ 
APE3 
HYP1 ! 
HYP3 
EXP1 

0-4 EXP3 
HEX1 > 
HEX3 



-0.2 -0.1 



-0.15 -0.1 -0.05 



0.05 0.1 0.15 0.2 



0.5 i 1 

APE1 I — f 
APE3 

HYP1 !-*- 
HYP3 i □ 
EXP1 
0-4 -EXP3 

HEX1 : ------ 

HEX3 



APE1 ^ 
APE3 
HYP1 i 
HYP3 
EXP1 
EXP3 
HEX1 i 
HEX3 



-0.1 -0.05 



0.05 0.1 0.15 0.2 



Figure 2: Our data for 2m PCAC vs. m with csw = (left) and Csw = l (right) at three couplings. 
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with m w = (l + b m am w )m w . The unmixed densities/currents X with X = S,P,V,A have 
been given in ffT2"j) . Note that either r.h.s. is scale-independent, since Z m — 1/Zg and the two 
renormalization factors Z$ and Zp run synchronously. Finally, due to the b m term in (iT4"|) . m cr it 
does not drop out of the r.h.s. of (TTTT) for unequal current quark masses. 

A naive determination of — am cr - lt = 4 — 1 / (2/c cr it) would measure M 2 as a function of tuq 
and determine, via an extrapolation, where the former vanishes. To avoid finite-volume and/or 
chiral log effects, we determine m PCAC as a function of uiq and see where this quantity vanishes. 
Up to 0(ag 2 , a 2 ) [depending on the details of improvement] cut-off effects this is a linear 
relationship, and, by virtue of (fT8l) . the slope is proportional to Z m ZpjZ^. More specifically, we 
restrict ourselves to degenerate quark masses (i.e. m 1 =m 2 ) and employ the fitting ansatz 

am PCAC = J- [1 + b m (am - am CTit )](am - am crit ) (19) 
Za 

with mo the bare fermion mass given in (TT3T) . The goal is to test how well the fitted — am cxlt 
and Za = ZA/{Z m Zp) agree with the 1-loop prediction. In principle, the coefficient b m is known 
at tree level. It turns out that using this value leads to unacceptable fits. On the other hand, 
our data are not precise enough to allow us to use b m as a parameter. The quoted fits use 
b m = 0; this leads in most cases to acceptable chisquares, and the few exceptions might be due 
to our limited statistics (cf. Tab. [7]). In fact, our data (taken at fixed am PCAC to limit the CPU 
requirements) do not show any visible curvature - Fig. [2] shows the data for three (out of five) 
couplings. We performed several alternative fits (e.g. by dropping the last data point), and as 
a result we estimate that the theoretical uncertainty is roughly one order of magnitude larger 
than the statistical error quoted in Tabs.lBlfTTl 

Our non-perturbative data for — am crit are given in Tab.[8]and Tab.[9]for the Wilson (csw = 0) 
and clover (csw = 1) case, respectively. As an illustration, we add the 1-loop prediction that 
follows from (0 [13]) and Tab.[TJ We did not measure the unfiltered — am crit , since it would be too 
expensive for our computational resources, and the large discrepancy between the perturbative 
and non-perturbative critical mass for unfiltered actions is well known. 

Our non-perturbative data for Za are given in Tab.fTUl and Tab. [Til for the cases csw = 
and csw = l, respectively. Note that Za is scale- independent, since Z m — 1/Zs, and the factors 
Zs and Zp run synchronously. Again, we add the 1-loop prediction that follows from ([7]) and 
Tabs.[5]-IU For similar reasons as above, we did not measure the unfiltered Za- 

The overall impression from Tabs. [MH] is that 1-loop perturbation theory does not give very 
accurate predictions for non-perturbatively determined renormalization factors, if the improve- 
ment coefficients are taken at tree-level. However, the mismatch is much smaller if filtering and 
improvement is used - as soon as one of the ingredients is missing, the "agreement" gets much 
worse. The virtue of the combined "filtering and improvement" program is that all renormal- 
ization factors and improvement coefficients are close to their respective tree-level values. This 
is in marked contrast to other schemes (e.g. [TB]) in which these quantities are far from and 1, 
respectively, and the challenge is to reproduce these big numbers in perturbation theory. 



4.2 Rational fits for m crit with APE/HYP/EXP/HEX filtering 

We know from ([5]) that asymptotically — am crit — * g^S/ (I2n 2 ) = S/ (27r 2 /3) with 5 given in Tab. CD 
Accordingly, if we fit our data with the rational ansatz 

- am crit = Cl f { ± C29 J (20) 

then the coefficient c\ would correspond, in the weak coupling regime, to 5'/(127r 2 ) with S given 
in Tab.[TJ Our data are not in the weak coupling regime, but still it is interesting to check how 
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c sw = 


P = 5.846 


P = 6.000 


P = 6.136 


P = 6.260 


P = 6.373 


thin link 


0.44573 


0.43429 


0.42466 


0.41625 


0.40887 


pert. 
1 APE 
lEXP 


0.11750 
0.5150(18) 
0.5846(23) 


0.11448 
0.4283(17) 
0.4932(15) 


0.11194 
0.3779(17) 
0.4412(21) 


0.10973 
0.3496(22) 
0.4039(11) 


0.10778 
0.3248(15) 
0.3805(09) 


pert. 
3 APE 
3 EXP 


0.04170 
0.2939(20) 
0.3263(21) 


0.04063 
0.2247(16) 
0.2509(18) 


0.03973 
0.1935(14) 
0.2100(19) 


0.03894 
0.1685(16) 
0.1869(08) 


0.03825 
0.1555(17) 
0.1713(19) 


pert. 
IHYP 
IHEX 


0.06046 
0.3094(18) 
0.3985(19) 


0.05891 
0.2455(14) 
0.3158(18) 


0.05760 
0.2093(16) 
0.2715(12) 


0.05646 
0.1908(15) 
0.2449(16) 


0.05546 
0.1728(28) 
0.2244(06) 


pert. 
3 HYP 
3 HEX 


0.1841(18) 
0.1993(18) 


0.1290(14) 
0.1419(14) 


0.1061(11) 
0.1142(16) 


0.0949(17) 
0.0976(14) 


0.0794(15) 
0.0868(16) 



Table 8: For Csw = Wilson fermions: — am crit with 8 filtering recipes. In each field, the first row 
gives the (common) 1-loop prediction, and the next two the linearly extrapolated values with 
APE/EXP or HYP/HEX filtering, respectively. Errors are statistical only. We did not measure 
the unfiltered — am crit , and we do not have a perturbative prediction for 3 HYP/HEX steps. 



c S w = l 


P = 5.846 


P = 6.000 


P = 6.136 


P = 6.260 


P = 6.373 


thin link 


0.27719 


0.27008 


0.26409 


0.25886 


0.25427 


pert. 
lAPE 
lEXP 


0.04254 
0.2438(13) 
0.3140(13) 


0.04145 
0.1929(08) 
0.2547(11) 


0.04053 
0.1685(07) 
0.2231(08) 


0.03973 
0.1518(08) 
0.2022(06) 


0.03902 
0.1413(05) 
0.1873(04) 


pert. 
3 APE 
3 EXP 


0.00668 
0.0779(13) 
0.1003(13) 


0.00651 
0.0497(07) 
0.0657(07) 


0.00637 
0.0400(04) 
0.0512(06) 


0.00624 
0.0341(02) 
0.0440(03) 


0.00613 
0.0312(03) 
0.0392(02) 


pert. 
IHYP 
IHEX 


0.01719 
0.0885(10) 
0.1464(14) 


0.01675 
0.0620(05) 
0.1045(07) 


0.01638 
0.0517(05) 
0.0851(04) 


0.01605 
0.0475(03) 
0.0743(03) 


0.01577 
0.0441(03) 
0.0674(04) 


pert. 
3 HYP 
3 HEX 


0.0252(12) 
0.0289(10) 


0.0120(05) 
0.0143(05) 


0.0094(02) 
0.0111(04) 


0.0084(01) 
0.0088(02) 


0.0077(03) 
0.0088(02) 



Table 9: For csw = l clover fermions: — am crit with 8 filtering recipes (cf. caption of Tab.EJ). 
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c S w = 


(3 = 5.846 


(3 = 6.000 


(3 = 6.136 


(3 = 6.260 


(3 = 6.373 


thin link 


0.94668 


0.94805 


0.94920 


0.95020 


0.95109 


pert. 
1 APE 
lEXP 


0.99859 
1.081(12) 
1.037(16) 


0.99863 
1.115(11) 
1.086(11) 


0.99866 
1.115(12) 
1.110(13) 


0.99868 
1.134(16) 
1.095(08) 


0.99871 
1.132(09) 
1.121(10) 


pert. 
3 APE 
3 EXP 


1.00281 
1.066(13) 
1.067(14) 


1.00274 
1.114(11) 
1.132(12) 


1.00268 
1.149(10) 
1.117(12) 


1.00262 
1.123(10) 
1.120(05) 


1.00258 
1.130(10) 
1.140(10) 


pert. 
IHYP 
IHEX 


1.00061 
1.046(12) 
1.070(12) 


1.00059 
1.129(09) 
1.117(12) 


1.00058 
1.120(10) 
1.126(08) 


1.00057 
1.137(09) 
1.139(09) 


1.00056 
1.114(14) 
1.129(07) 


pert. 
3 HYP 
3 HEX 


1.051(11) 
1.058(11) 


1.113(09) 
1.119(09) 


1.119(08) 
1.125(10) 


1.148(11) 
1.125(10) 


1.114(10) 
1.119(09) 



Table 10: For Csw=0 Wilson fermions: Za/ (Z m Zp) with 8 filtering recipes (cf. caption of Tab.EJ). 



c S w = l 


(3 = 5.846 


(3 = 6.000 


(3 = 6.136 


(3 = 6.260 


(3 = 6.373 


thin link 


0.90710 


0.90949 


0.91149 


0.91324 


0.91478 


pert. 
1APE 
lEXP 


0.98030 
0.8523(86) 
0.8602(80) 


0.98080 
0.9263(55) 
0.9477(44) 


0.98123 
0.9679(51) 
0.9912(26) 


0.98160 
0.9798(70) 
1.0079(17) 


0.98193 
0.9974(42) 
1.0200(24) 


pert. 
3 APE 
3 EXP 


0.99424 
0.8554(69) 
0.8458(80) 


0.99439 
0.9398(33) 
0.9503(28) 


0.99451 
0.9761(31) 
1.0024(15) 


0.99462 
1.0008(22) 
1.0258(10) 


0.99472 
1.0126(24) 
1.0340(21) 


pert. 
IHYP 
IHEX 


0.99014 
0.8539(88) 
0.8703(78) 


0.99040 
0.9195(72) 
0.9510(42) 


0.99061 
0.9599(51) 
0.9841(34) 


0.99080 
0.9699(49) 
1.0039(25) 


0.99096 
0.9847(37) 
1.0160(15) 


pert. 
3 HYP 
3 HEX 


0.8517(93) 
0.8452(61) 


0.9435(42) 
0.9548(30) 


0.9713(33) 
1.0050(26) 


0.9959(23) 
1.0193(13) 


1.0074(28) 
1.0373(11) 



Table 11: For csw = 1 clover fermions: Za/ {Z m Zp) with 8 filtering recipes (cf. caption of Tab. El). 




APE1 

APE3 — 
HYP1 x 
HYP3 : n 
EXP1 
EXP3 
HEX1 • 
-HEX3 



' 

0.94 



Figure 3: — am cri t vs. g% for Wilson (csw = 0, left) and clover (csw 
filterings. The curves indicate 3-parameter rational fits. 



1, right) fermions with 8 
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c sw = 


csw = 1 


pert. 


n 1 1 A a on 


U.U41440 ( 


1 APE 1 EXP 


0.213(12) 0.252(12) 


0.0909(28) 0.1094(20) 


pert. 


U.U4Ubzy 


U.UUooUyo 


3 APE 3 EXP 


r\ r\TT/~\ a\ r\ aoo / r\'~7\ 

0.077(14) 0.083(07) 


0.0172(15) 0.0171(09) 


pert. 


U.UooyUo 


U.UlO f OUz 


1 HYP 1 HEX 


0.095(14) 0.121(04) 


0.0338(12) 0.0332(16) 


pert. 
3 HYP 3 HEX 


0.034(15) 0.026(01) 


0.0060(02) 0.0060(15) 



Table 12: The fitted coefficient c\ in (1201 . compared with the 1-loop prediction S'/(127r 2 ) with 
5* taken from Tab.[TJ 



APE1 

APE3 - 
HYP1 i * 
HYP3 n 
- EXP1 
EXP3 

HEX1 !.-..- 
•li X3 



0.94 0.95 0.96 0.97 0.98 0.£ 




1.01 1.02 1.03 




1 1.01 1.02 1.03 



9o 9o 

Figure 4: Z A /(Z m Z P ) vs. g 2 for Wilson (csw = 0, left) and clover (csw = l> right) fermions with 
8 filterings. The curves indicate 3-parameter rational fits. 



much the coefficient c\ from an unconstrained fit deviates from the perturbative value. The 
result is shown in Fig.[3]and Tab.fT2l As was to be anticipated from our discussion of TabsJEHHJ 
the "agreement" is not very good. On an absolute scale the numbers are close, since they are 
all much smaller than one. On a relative scale, they deviate by a substantial factor. In spite of 
this disagreement, the non-perturbative data still show a consistency c^ PE ~cf xp and ditto for 
c^ YP ~c™ x , as predicted in PT. We find this amusing, in particular in view of the fact that the 
corresponding (in PT) curves in Fig. [3] are not close at all. 

4.3 Rational fits for Z A with APE/HYP/EXP/HEX filtering 

We know from ([7]) that asymptotically Za — > 1— g^ZA+zs— zp)/{12'K 2 ) = l—{zA+zs—zp)/{2n 2 (3). 
Accordingly, if we fit our data with the rational ansatz 

~ 1 + d ig 2 + dggp . . 

1 + d 3 g 2 

then d\— d 3 would correspond, in the weak coupling regime, to (za+zs~ zp)/{12ii 2 ) with za, z$, zp 
given in Tabs.[2]-Hl The result of our fits is displayed in Fig.Hl Again, there is no quantitative 
agreement between 1-loop perturbation theory for Za and our non-perturbative data, based on 
tree-level improvement coefficients. Still, comparing the two graphs in Fig.Hl one is l e d to believe 
that with appropriate 1-loop improvement coefficients the situation might be better. 
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4.4 Rational fits for m res with APE/HYP/EXP/HEX filtering 

We may express our result in terms of m res = m PCAC (m = 0). Z& ~ 1 implies m res ~ — m CT i t , 
and we refrain from copying Tabs. [BE with minimal modifications. Again, we performed ra- 
tional fits, and the result looks very similar to Fig. [31 An interesting observation is that 
m res in physical units is almost constant. We find m^ PE ~ 144, 111, 107, 108, 113 MeV at 
(3 = 5.846,6.0,6.136,6.260,6.373 and m™ YP ~ 47, 27, 25, 26, 27 MeV. We feel confident that 
with 1-loop values for the coefficients csw> c a, b^—bp smaller residual masses could be obtained. 



5 Summary 

We have presented a systematic study of filtered Wilson and clover quarks in quenched QCD. 
We have derived results at 1-loop order in weak-coupling perturbation theory for — am crit and 
the renormalization factors Z x with X = S,P,V,A with four filterings [APE, HYP, EXP, HEX], 
in some cases with 1,2,3 iterations. We have compared these predictions to non-perturbative 
data for —am cr \ t and Za = ZaZs/Zp in a simulation without improvement and with tree- level 
improvement coefficients. We find no quantitative agreement in this specific setup. Still, the 
tremendous progress that comes through the combination of tree-level improvement and filtering 
leaves us optimistic that a theory with 1-loop improvement coefficients and 2-loop renormaliza- 
tion factors might work in practice. By this we mean that a continuum extrapolation can be 
done from accessible couplings as if the theory would have 0(a 2 ) cut-off effects only. 

It turns out that lattice perturbation theory for UV-filtered fermion actions is not much more 
complicated than for unfiltered actions. For instance, our formula f)73p gives a compact 1-loop 
expression for the critical mass with an arbitrary number of APE smearings, and shows that 
Gw^crit - ► for njtor oo. Since our results in the main part of the article were derived in a fully 
automated manner, we feel that this explicit calculation provides an important check. 

One particularly compelling feature of filtered clover actions is that tadpole resummation 
is not needed; in fact it barely changes the result. This suggests that perturbation theory for 
filtered clover quarks converges well. In consequence, we expect that for filtered clover fermions 
the non-perturbative improvement conditions as implemented by the ALPHA collaboration [15] 
will yield values consistent with such perturbative predictions. 

A beneficial feature in phenomenological applications is the low noise in observables built 
from filtered clover quarks. We have been able to determine m crit to ~ 3% statistical accuracy 
from just a handful of configurations. Therefore, the "filtering" comes at no cost - it actually 
reduces the CPU time needed to obtain a predefined accuracy in the continuum limit. 

Let us comment on the filtering in two different fermion formulations. It is clear that twisted- 
mass Wilson fermions would benefit from filtering, too. The dramatic renormalization of the 
twist angle would be tamed and it would be much easier to realize maximum (renormalized) twist. 
For rather different technical reasons, filtering has proven useful for overlap fermions [IHIES> 30J. 
In our technical study we decided to stay with csw = 0, because the overlap prescription achieves 
automatic 0(a) on-shell improvement. It is not clear to us whether the better chiral properties 
of a clover kernel could translate into further savings in the overlap construction. 

We hope that, once the 1-loop value for csw with n iterations of the EXP/stout recipe [20] 
is known_|, filtered clover fermions are ready for use in large-scale dynamical simulations. An 
important point is, of course, the smallest valence quark mass that can be reached for a given 
coupling and sea quark mass (partially quenched setup). We find am^^ YP — 0.0126(5) at [3 = 6.0 
and am r i e ^ YP = 0.0074(3) at /3 = 6.373 in the quenched theory. This corresponds to an almost 

3 Note that at 2-loop order the strict correspondence between APE and EXP with a APE /6 = a EXP is lost. 
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constant residual mass in physical units, mj? e ^ YP ~27MeV. Since this mass is much smaller than 
in the unfiltered case, it is natural to hope that one can reach smaller valence quark masses (in 
the quenched or partially quenched setup) before one runs into the problem of "exceptional" 
configurations. Furthermore, if mixing with unwanted chiralities in 4-fermi operators is an 
0((am res ) 2 ) effect [31] in our case, too, the small residual mass would be relevant for electroweak 
phenomenology. Clearly, these topics deserve detailed investigations. 
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A Fat link perturbation theory in d dimensions 
A.l APE smearing 

In d dimensions and with general gauge group G, standard APE smearing is defined through 
U' fl (x) = P G {(l-a)U^x) + -^— £ U p {x)U,{x+0)U p {x+^\ 

where the sum ("staple") includes 2{d— 1) terms. The projection P G is needed, since in general 
the staple is no longer a group element. For the perturbative expansion we substitute U^(x) — ► 
1 + iaA fl (x+ 1) + 0(a 2 ). The prefactors 1 — a, a/(2d— 2) ensure that in PT the effect of P G is 
already taken care of. For 2-quark and 4-quark renormalization factors at 1-loop order only the 
linear part is relevant [11]. After shifting x — > x — | one obtains at leading order 

a 



A'^x) = A^x) + — —^{A^x+0)-2A ll (x)+A lt (x-0)} 



+ -r, :V \A V + - -AJx-- — - AJx+-+-) + AJx+- — 1 22 

2(d-l) A ^^ v 2 2 } y 2 2 1 v 2 2 } v 2 2 J i v ' 



2(d 

where the sum now extends over all positive v. This may be recast into the form 



Vf+f Vf-| + 5 v 



a 

a 

2{d-l) 

A'^x) = E KM A„{x+y) = J2 {Ks/)<W + «V(l/)] A v (x+y)} (23) 
which is suitable for a Fourier transformation. This leads to the final relation 

at . . . .. , . . a 
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with q p — | sin(|g p ) (for all d). A form particularly useful for iterated smearing {n> 1) is [TT] 

4%) = E {([i - ^iy^ n (<W - ^r) + ^r) M*)} (25) 

where the transverse part is simply a form-factor f^(q 2 ) = [1 — 2( -Jl 1 - ) g 2 ] w as emphasized in 
A. 2 HYP smearing 

In d > 3 dimensions d— 1 levels of restricted APE smearings may be nested in such a way that the 
final "fat" link contains only "thin" links in the adjacent hypercubes [9]. Specifically, in d — A 
the linearized HYP relation reads (note that 0:3,2,1 refer to step 1,2,3, respectively) 

A^ P (x) = (l-a 3 )4(^) + y £ {A a {x-^ + ^)+A fl {x+a)+A cr {x+^ + ^)} 
A^(x) = (l-a 2 )A lx (x) + ^ J2 { A P^(x-^ + ^)+A lx>vp (x+p)+A Pjlll/ (x+^ + ^)} 

Afa) = (l-ai)4(x) + y {AuA x -% + ^)+AvA x +v)+A u A*+% + %)} 

and it is easy to see that the core recipe in each step is an APE smearing in 2,3,4 dimensions, 
respectively Therefore, the Fourier transform leads to the relations 

A^ P (q) = 4(<z) + y E {-fiM<i) + w*MQ)} 

A^(l) = (1 - "2)4,(9) + E {(2-9p)4^(?) + ^4,^(9)} 
4(g) = (l-ai)^( g ) + ^{(2-ai,(?)+?M,^)} ( 26 ) 

where a simplification specific to the innermost level has been applied. Plugging everything in 
we obtain a compact momentum space representation for one level of HYP smearing 

4 = 4 + y E {»A,-<?'4 + ^ E {29^4- ^[(2-^)4 + ^^] + 

^ ]T {4g^4.-g'[2g M gpA p + (2-®[(2-^ + g M M.]]}}} (27) 

which, however, entails some orthogonality constraints. To get rid of the latter, we apply a 
number of tricks. First, the sum over o is split into two parts. The part quadratic in q a can 
be made independent of the summation index by virtue of q 2 = q 2 — q 2 , — q 2 . — q 2 .- Hence, what 
remains in the innermost summation is the term linear in q a . This term, however, is independent 
of p, the next-level index. Since the constraint lets it assume the other free value (after p and v 
have been fixed) than p, the total effect is the same as with a — > p replaced, thus 

4 = i + + 7 E {(2+«3(2-g 2 + gJ + ^ + $)g^4 

+ f (f -ft- ft- ©(2 - <©H(2 - €4 + UM}} 
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is a representation with only two sums. Next we pull out those parts which are independent of 
the index p. Using Y.p^v) g 2 — <? 2 ~~ Qu ~ ql m the remainder yields 

4 = + y E + x E {( 2 +«3(2-g 2 + ^ + C + ©)»4} 

-f [(i+« 3 )(? 2 -c-a-y[(? 2 -c-^) 2 - E 5»2-a^+wA]} 

where the bracket multiplying ^ is just 2 ]lp^(pV) g 2 - Since a constrained product would be 
inconvenient for later use, we choose to stay with the actual form, but now we relax the constraint 
on p to differ from p only and compensate for the additional term. This yields 

4 = ^m + y E {[l-f (2 + « 3 (2-g 2 + gJ + 2a)]»^-?X 

"f [(1 + «s)(g 2 - Ql - €) ~ y«(2 - ql)A, + %q v A v \ 
+^r£ E KW + f £ {(2+«3(2-g 2 + gJ + g 2 ))g^}} 

with Qpj, = (g 2 — q 2 — ql) 2 — Z)p^(^) <Zp- I n the sum over p the term which depends on v has 
been isolated. The reason is that the other term may be pulled out of the z/-sum (this yields a 
factor 3), and since the constraint is the same, renaming the index p^>v is then legal. Applying 
a similar procedure to the ^-independent factor of the former term, we obtain the form 

4 = 4 + y E {[i + Mi + j(4-f + g; + C)E?A-?X 



—la + «3)(<f - c - a - f«(2 - <? 2 m, + 



with just one sum [apart from the A-independent J2 q A p in Qy,v = (g 2 ~qfi~ ql) 2 + <?p + qt ~ Sp Qp] 
Now it takes a couple of algebraic manipulations to arrive at the form 

4 = 4 + y E O 1 + "2(1 + « 3 - |(i + 2« 3 )(g 2 - g 2 - g 2 ) + - g 2 ^ 

[(1 + a 3 )(2<f - 2gJ - 2q 2 u ) - a 3 ((f) 2 - 2q% + 2gJ - £ g p 4 ) - (1 - a s ) x 

(g 2 - ql - €)€ + f((g 2 ) 2 - 2g 2 gJ - mi + K + + 2 <?' - E 

which is suitable to do the sum in the terms which are even in q v . This operation yields 

4 = 4 + y E {[i + «2(i + a3-^(i + 2« 3 )(g 2 -gJ-g 2 ) + yQ^)]g^g^} 

- y(l + « 2 (1 + « 3 ))(g 2 - ©4 + ^y (1 + 2« 3 )((g 2 ) 2 - 2g 2 gJ + 2gJ - £ #4, 
aiQ;2Q; 3 



(((g 2 ) 2 - 2<f <g + 4gJ - 3 E g p 4 )(g 2 - g 2 ) + 2 £ g p 6 - 2gJ)A, 



48 

and upon extending the sum and compensating for the additional term one finds 

4 = ^ + yE{[ 1 + « 2 (l + « 3 -^(l + 2« 3 )(g 2 -g;-g 2 ) + f Q,u)}%q»A v } 
-y [1 + «2(1 + « 3 )]g% + ^(1 + 2a 3 )[(g 2 ) 2 - g 2 ^ - £ 

ai«20 3 r/~2\3 



- 2(g 2 ) 2 gJ + 2g 2 ^ - 3g 2 £ g p 4 + 2gJ £ g p 4 + 2 £ $4 (28) 



48 LVJ 
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which looks somewhat lengthy. As was noted by DeGrand and collaborators [TTJ [TE1 EH] , defining 
Qfiv = 1 + 012(1 + a 3 ) — ^-((1 + 2a 3 )(g 2 — q 2 — q^) — ^-Q^ v ) allows for the compact form 

4 = E {(! - f ERpOH* + f (29) 



without any constraint on v or p. The general form for iterated smearing (n>l) is 



4° = £{r£?(* 



Q 



(30) 



with the transverse and the longitudinal form-factor both being the product of n factors with 
adjacent indices summed over and the first and last index set to \x and v respectively, 



t(") 
p^ 



L (n) 



e n(i-|a[«^+«rf«} 

Ai,...,A n _i i=l Pi 



A =p,A n =^ 



e n (1 -la^+^iai + ^jv^ 



6 



Ao=/i,A n =^ 



(31) 
(32) 



Ai,...,A n _i i=l " Pi 

In practice only moderate n are relevant, and for n = 2 and n = 3 the explicit formulae read 



4 2) 



E{(i-yE{«)V + 



6 „ 

,ll ^( 2 - f E{[^P+^p]?p 2 }) + ^E{^p^?})^}^ 



6 



4 3) 



E{(i-fE{^})V + (> 



Oil 



3-^E{[^p + ^p]?p 2 } 
z p 



+ |t(E{[^p + ^p]<?p 2 }) 2 - ^E{^p 2 >E{^a€} 



36 



2^ E^P^^A^a})^^ 

P,A 



(33) 



(34) 



but it is still clear that in general the transverse part contains a factor (1 — ^ Ep{^pp<?«}) n 



A. 3 EXP smearing 

Here we consider the EXP/stout smearing U'{x) = S fl (x)U fl (x) [no sum] introduced in [20; with 

S,(x)=exp(^{[ £ U v (x)U^x+u)Vt(x+fi)Ui(x)-U^x)U v (x+^ 

A ±^(p) 6 

a special unitary matrix by construction. Upon expanding as before we obtain 
l + iaA^a;) = (l + iaa £ {A v (x-^+^)+A^(x+u)-A v (x+^+^)-A lt (x)}){l + iaA li (x)) 



±^(p) 
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and thus (still, up to terms of order 0(a ,)) 

A'^x) = {l-2(d-l)a)A li {x) + a £ + A^x+fi) - A u {x+^+^)\ 

±u^{p) A A / A 

= A„(x) + aJ2{A f ,(x+i))-2A f ,(x) + A^x-i>)} 



+ {M*-% +£) - M^~~) - M*+%+%) + M^~)} (35) 

which is just (1221) with a modified parameter. Accordingly, 1-loop fat link perturbation theory 
for EXP/stout smearing follows from the version for APE smearing through the replacement 

* 2(d-l)a EXP / stout . (36) 



a 



APE 



A. 4 HEX smearing 

A natural generalization of the HYP concept is to use EXP / stout smearing in each of the 3 steps 
(in 4D) rather than the standard APE smearing [SJ. This entails the general definition 



V,, vp {x) =exp(^{[ £ l^ n - 1 H*)^ 1 H*+^)^ 1 H*+A) t ^ ,1 - 1) (*) t - h.c] - 

±<r£(p,v,p) 

%, v {x) = exp(^{[ £ ^(^)^,-(^ + ^)K^(^+A) t ^r 1) (^) t -h.c.] --&[•]}) 
U^\x) 

U { ;\x) =exp(^{[ £ V^(x)V^(x+u)V^(x + 



2 ±v*(p) 



U^\x) (37) 

where again 03,2,1 refer to step 1,2,3, respectively, and no summation over \x is implied. We refer 
to (EZD as "hypercubically nested EXP" or "HEX" smearing. With ([36]) it follows that 

^HYP „HYP „HYP\ (fir. HEX a „ HEX r> HEX\ / Q o\ 

(a 1 ,a 2 ,a 3 ) — > (6a 1 , 4a 2 , 2a 3 ) (38) 
will automatically generate the perturbative formulae for the HEX recipe (!37l) . 

A. 5 Permissible parameter ranges 

Regarding a reasonable range of smearing parameters, a standard criterion that one may impose 
to avoid instabilities at higher iteration levels is that the form-factor shall be smaller than 1 in 
absolute magnitude over the entire Brillouin zone. Since q 2 < 4d, formula ( |25l) gives 



a max j (39) 



for APE smearing with arbitrary iteration number n. With the replacement prescription (1361) 
the analogous condition for EXP/stout smearing is a EXP < 5-3. 

For n HYP smearings in 4D the transverse part contains the factor (1 — ^ Ylp{^ppQp}) n i an d 
requiring this to be bounded in absolute magnitude by 1 leads to the two-fold condition 

< Z{ ai (l+a*(l+a*)-%[(l+2a«)(?-^^ < 12 

P 4 Z A 
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for each p. Accordingly, upon summing everything over p one finds 

< ^{ ai (4 + 4« 2 (l + « 3 )-f [(l + 2a 3 )(3g 2 -4g 2 )-^ [2(<f ) 2 -6«f 9 2 + 8^ -2 £ gl]])g 2 } < 48 

P 4 Z A 

and then doing the sum over p yields the inequality 
0<4e* 1 (l + a 2 (l + a3 ))<f-^[(l + 2^^ 

4 A A A 

which is a non-trivial constraint on (ct^ YP , a^ YP , ) in terms of the three quantities 
0<£<?a<i6 , 0<^gl<64 , 0<£$<256 

AAA 

but the latter are, of course, not independent. Neglecting this, the condition 

< ai(l + a 2 (l +a 3 ))[0...64] + «ia 2 (l + 2a 3 ) [-192...64] + aia 2 a 3 [-1024...1280] < 48 (40) 

can be separated into one for the lower and one for the upper bound. While the former is always 
satisfied for positive smearing parameters, the latter takes the form 

ar P (l + a» YP (2 + 23af YP ))<^. (41) 

Another useful form might arise from keeping only the part quadratic in the momenta in the 
inequality, as the remainder may have either sign, and this leads to the less restrictive condition 

HYP/i | „HYPh i ^HYP\\ ^ 3 . . 

a x [l + a 2 [l + a 3 )) < - . (42) 

Note that for a 2 = a 3 = either condition coincides with (139]) . Finally, we mention that neither 
(jHj) nor (|42|) is satisfied by the standard HYP parameter set PJ. Note, however, that these are 
not necessary conditions; they emerged from applying some simplifications to a highly non-linear 
precessor. Applying the replacement recipe (1381) . the analogous conditions for HEX smearing are 
found to be a™ x (l + 8a™ x (l + 23a™ x )) < § and af EX (l + 4a 2 IEX (l + 2a™ x )) < |, respectively. 

A. 6 Diffusion law for iterated smearing 

As a consequence of (1251) . the form-factor for the transverse part after n APE smearings is [11] 

f n \f) ~ ex P (-^ ri yg 2 ) + 0((g 2 ) 2 ) . (43) 

This means that the square-radius of the resulting form-factor takes the form 

na APE 

(Oape = (44) 

which is a diffusion law, since the smearing effectively affects a space-time region growing like 
( r2 )APE 0C V / ^'- Focusing on the quadratic part in the transverse factor below fl34l) one finds 

(r 2 >HYP = ^ i (l + a 2 (l + a 3 )) (45) 



for n iterations of HYP smearing in 4D. As noted in [11] . the prefactors are favorably small. 
Even 3 APE steps with a^ PE generate a "footprint" (?" 2 )ape — 0.775 i.e. of the order of one 
lattice spacing. Likewise, 3 HYP smearings with a^ YP yield (r 2 )JjYP — 1-155. 
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B Additive mass renormalization with filtering 

Here we give a derivation of the additive mass renormalization for APE-filtered clover fermions 
at 1-loop order in lattice perturbation theory. We work in Feynman gauge; the effect of smearing 
is not just a modification of the gluon propagator, as it is in Landau gauge [TT] . 
For the gauge field we use the same conventions as in App. A, that is 

4%) = h${q)Af{q) (46) 

= (l^)1^-% + ^ = /m,-(f(9)-# (47) 
6 q A q z q z 

with f(q) = 1 — (a/6)q 2 and (7 = 2sin(<7 At /2), except that repeated indices are always summed over 
in this appendix. Furthermore, we use the shorthand notation 

. / q^ \ 2 

= sm(q p ) , s 2 = s^s^ 

and analogously c At = cos(g /i /2) and c 2 = c At c At with summation implicit. 

With these conventions the gluon and quark propagators (in Feynman gauge) take the form 

G, u (q) = 5, u G{q), G( ff ) = ^ (48) 

S{q) ~ A(q) - A(s 2 ) 2 + P (49) 
and the two-quark (zero external momentum on one side) one-gluon coupling is V P ±W P with 

V P (q) = -ij P Cp-s p (50) 
W p (q) = -^r-(yp\c p sx (sum over A only) (51) 

where we have separated the csw independent part from the part linear in the clover coefficient. 
The precise form of ( 1501 15TI) refers to the U(l) gauge theory; we will include a factor Cp below. 

B.l Sunset diagram 

With V} n ) = h p ^ V a the part of the sunset diagram proportional to (ap)° follows from 

ls ^ t]MCF) = jj^ Giq) m<)+wm<)-w?>\ (52) 

v-mbv-w = fxBK - nr - i) s > 8 ° v ° BV > ± "j^sMYs. + (r _ i f>" s ^f Vi > 

s [s ) 

= f-VpBV p + {l-f 2n ) S -^V a BV p (53) 

and analogously for V^BW^ n \ W p {n) BV p {n) and W^BW^l The terms even in q are 

V a BVp = 2s 2 s a S/3 - 2^ a ^f 3 c a Cf } s 2 + (7„7 A1 c a s /3 + ^^ps a cp)s p (54) 
V a BWp = —(i a 'y p af]\c a c (3 s IM sx + 2a l3 xs 2 cps a sx) (55) 



2i 

W a BV/3 = ( ^-(a a xl^lpc a cps Pj s x + 2a a xs 2 c a Sf 3 s x ) (56) 



c 2 

W a BW p = — ^cr^a^xCaC^s^x (57) 
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where = stands for "up to terms odd in q" . With this at hand, we compute 



V P BV P 
V P BW P 



W P BV P 



W P BW P 



s a s p V a BVp 
Sf3 V a BW p 
s p W a BVp 

80 W a BW p 



2(s 2 ) 2 -2(A-s 2 )s 2 + s 2 

{l P l^p\C 2 p s p s x + cr pX s 2 s p s x ) = first + 



CSW/ 2 ■>- - 

( " ,.a. \c ■■ ■ - 

2i 

csw 
2i 



(vpxlfilpC p s p s x + (r pX s s p s x ) = first + 



--SW 2 2-- 

— OpnOpXCpS S K S X 



2( S 2 )3 + V-2 

^-(jlal^pXSaSpS^Sx + O p X (s 2 ) 2 's pS X ) = 

C SW/1 . / 2\2 \ n 

-^-{^aXlulpSaSpS^Sx + (T aX (s ) S a S X ) = 



-sw 



v aK vpxs S a SpS K S X = 



(58) 
(59) 

(60) 
(61) 
(62) 
(63) 
(64) 
(65) 



where the asserted vanishing of certain terms holds only in case there are no further factors 
which destroy the symmetry property it builds on. With fl53l we thus arrive at 



BV [n) 



V p {n) BWl n) 



W^BW^ 



f 2n {A(s 2 ) 2 - 8s 2 + s 2 ) + (1 - / 2n )(2(s 2 ) 2 + ^s 2 ) 
1 



2(8 



2\2 



^-s 2 + f 2n (2(s 2 ) 2 - 8s 2 + ~s 



7:2 \ 



2i 



W r-2n 2- - 

-/ lplp.VpXC p S p S X 



2i 



-/ n (y P xii,i P c p s p s x 



-SW r2n 



fin 2 2-- 

/ OpK^pXCpS S K S X 



(66) 
(67) 
(68) 

(69) 



making the numerator in ( |52l) take the form 



\y{n) + w^ n) }B[V p {n) - W^} 



+ 



+ 



2 ( S 2 ) 2 + ^2 + /2 n (2(s 2 )2 _ g s 2 + ± -2) 



-SW j2n 



f Zn v pK a pX c 2 p s 2 s K s x . 



(70) 



By means of the identities 0^7^ - jpj^a^, = 2i[y x -f p - 5 Xp 5 Pll ] and o pK o pX = j K -f X - j K j p 5 pX - 
bnplplx + finpfipx, where in either case p is not yet summed over, we thus obtain 



[sunset] /{glC F ) 



d 4 q f 2n - 



A f 



2n 



(2vr) 4 4(s 2 ) 2 + s 2 



+ c S w 



A f 2n 44 - c 2 p s 2 p 



,2\2 1 5 2 



+ 



-sw 



(2vr) 4 4s 2 4(s 2 r + 
A / 2 ^ 2 [c>l-c^] 



2 J (2vr) 4 4s 2 4(s 



2^2 1 q2 



where Z = Jd 4 q/(27r) 4 l/(4s 2 ) = 0.15493339... has been used. 



(71) 
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B.2 Tadpole diagram 

The tadpole diagram is readily evaluated to give 

d A q G{q) 



[tadpoleJo/OoCV) 



-4 



(2tt) 4 2 

'4, 



ha 



(with p fixed) 



1 f^Q_l T ( f n 5 

2 7 (2tt) 4 s 2 ^ V J °p« 



1 



(27T) 

d 4 g 



(r 



2 7 (2vr) 4 s 2 
1 r d A q 1 



/ 2n_ 2r(r _ 1) | + (r _ p2 9p 



q 



2 y (2vr) 4 s 2 

i 



r + a - / 



2n\ 





/■ A 3/ 2n i 


Z + 


7 (2tt) 4 4s 2 -1 



where Z = /dV(2vr) 4 l/(4s 2 ) = 0.15493339... has been used. 



(72) 



B.3 Combining the two 

It is now straightforward to add (1711) and (1721) to obtain for am cr \t = So the result 

d A q f 2n 



d*q f 2n 

+ 2 



(2tt) 4 4s 2 J (2vr) 4 4(s 2 ) 2 + s 2 
csw 



(2tt) 4 4s 2 4(s 2 ) 2 + s 2 



a 



sw 



d\ 



f 



r 2 ? 2 
2n C P S X 



c 2 s 2 



(73) 



(2tt) 4 j 4(s 2 ) 2 + s 2 

and we comment on the four contributions. The first term without the 1 / (4s 2 ) factor would be 

^ dki...dkd 



j(m) 



a 



2(d-l 



-k 



n (27l) d 
d m [*dki...dk d a[l-^-r 

j g 1 2{d— 1 

a= J-n (27T) d 



da 

d" 



da r 



cr=0 



era / (7 (X 

e d-i/n 



d- 1 



(74) 



with m = 2n, Iq denoting a Bessel function of the second kind and i^ -* = 1. The first term with 
the denominator but without the smearing factor would assume the simple form 



J(0) 



w dki...dk d 1 
(2vr) d k 2 



dr 



-2t 



/o(2r) 



f* dki...dk d & 
' w (2tt)^ 

oo (d = 2) 

0.25273101... (d = 3) 
0.15493339... (d = 4) 



(75) 



and the actual first contribution can hence be handled via a recursion formula 



j(2n) 



w dki...dk a 

-7T (2yr) d 



1 - 



a 



2(d-l) 
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2n 1 





n = n — 1 n = 2 n = 3 


C SW 
A 

c sw 

r 2 


51.43471 13.55850 7.18428 4.81189 

1 Q 70010 a Afil qo A 7r\A£,7 Q KfiflRC; 
lO.lOOio O.yoloo 4. rU40( O.OOUOO 

45.72111 13.50679 6.52280 3.84215 



Table 13: Numerical values of the integrals in (!73|) for a APE = 0.6 and n = 0..3 iterations. 



n dki...dk a 

j(2n-l) 



1 - 

a 



a 



2(d 



2(d-r 

.j(2n-l) 



■k 



2n-l 



1 



k 2 2{d-l) 



J(0) 



2(d 



j(0) + j(l) + + j(2n-l) 



(76) 



and ditto for 2n — > m. For the other terms we resort to numerical integration. We collect the 
pertinent values in Tab. [TBI With these it is easy to verify the APE entries in Tab.[TJ 



B.4 Other smearing strategies 

In this article we have focused on a strategy where one applies the same smearing in three 
places: in the covariant derivative and the Wilson term of the Wilson operator ([lj and in the 
field-strength tensor of the clover term (T5]). Of course, other options are possible. In general one 
may apply n steps with parameter a to build the links for the (relevant) covariant derivative, n! 
steps with parameter a' in the Wilson term and n" steps with parameter a" for the clover term. 
The numerator in ([52]) then takes the form [V^ n '\q) + W^ n '' ) }B(q)\V^ n '\q)-W^'' ) ] where n 
denotes the smearing level in the (relevant) covariant derivative, n' that in the Wilson term and 
n" the one in the clover term. Possible choices include: 

• n = n! = n" = 0: standard (thin-link) clover action (SC) 

• n = 0,n' = n" > 0: fat-link irrelevant clover action (FLIC), Wilson and clover terms [T3] 

• n > 0, n' = n" = 0: fat-link relevant clover action (FLRC), only covariant derivative 

• n = n' = n" > 0: fat-link overall clover action (FLOC), same smearing everywhere [T2"l[T4"] 

All explicit numbers given in this article refer to the "FLOC" case, but it is straightforward 
to generalize the formulae to arbitrary n, n', n". For instance, for n = n' the terms in (1731) 
proportional to c$ w , Cg W , Cg W contain a factor f 2n , f n+n , f 2n , respectively. 



C Details of the parameter dependence 

In this article we focus on the "standard" parameters ([3]) for APE/EXP smearing and (jlj) for 
HYP/HEX smearing. Here, we briefly discuss the dependence on a APE = 6a EXP . 

In Tab. fT4lfT6l we give details on how 5* and zx for X = S, P, V, A depend on the smearing 
parameter with 1,2,3 steps of APE/EXP filtering with csw = 1- I n most cases, one finds a 
reduction of S and (zp — z$)/2 = zy — ZA for a APE between and ~ 0.75; beyond that they 
increase sharply. This is in line with the discussion in App. B - perturbatively, one expects 
larger smearing parameters to be more efficient, up to 0^^ = 0.75 or a EXP = 0.125. Hence our 
"standard" choice ([3]) for the smearing parameter is not bad - at least in perturbation theory. 
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1 A T~»TH 

1 APE 


0.12 


0.24 


0.36 


r\ AO 

0.48 


0.6 


r\ to. 

0.72 


0.84 


0.96 


S 


23.51856 


16.57684 


11 1 /"» 1 O 1 

H.16131 


7.27195 


4.90876 


a /~\ m 1 T f 

4.07175 


4.76092 


6.97626 


zs 


15.01627 


11 O A f\ f A 

11.34954 


o o r\r\!~7 r* 

8.30976 


f o c\ r* c\ a 

5.89694 


4.H106 


2.95213 


2.42016 


2.51513 


Zp 


17.42350 


13.18607 


r\ /ifr r\r\ r\ 

9.67028 


6.87614 


a r\ r\ o r* a 

4.80364 


3.45280 


2.82360 


2.91606 


Zy 


11. 1 48/0 


o. 1 o69o 


6.3o36z 


4.o3o (a 


O 1 1 A O 

o.31z4o 


Z.b /4o6 


Z.DZOlo 


o 1 a a on 
0.164/y 


ZA 


10.54515 


7.83869 


5.67337 


4.04918 


2.96614 


2.42423 


2.42346 


2.96382 


Table 14: S and 


zx versus smearing parameter for 1APE clover fermions with csw = l- 


2 APE 


0.12 


0.24 


0.36 


0.48 


0.6 


0.72 


0.84 


0.96 


S 


17.57370 


9.36111 


4.99413 


2.77532 


1.66435 


1.27800 


1.89014 


4.43178 


zs 


11.77989 


6.90750 


3.79645 


1.77701 


0.40606 


-0.53293 


-1.02988 


-0.84811 


zp 


13.68411 


8.05861 


4.48219 


2.185H 


0.65185 


-0.37897 


-0.91453 


-0.70780 


zv 


9.15359 


5.43316 


3.29869 


2.10671 


1.43934 


1.10432 


1.13496 


1.79020 


ZA 


8.20148 


4.85761 


2.95582 


l. 90266 


1.31645 


1.02734 


1.07729 


1.72005 


Table 15: S and 


Zx versus smearing parameter for 2 APE clover fermions with csw — 1- 


3 APE 


0.12 


0.24 


0.36 


0.48 


0.6 


0.72 


0.84 


0.96 


S 


13.33394 


5.67072 


2.60949 


1.34010 


0.77096 


0.57512 


1.14061 


4.42535 


zs 


9.29121 


4.16584 


l. 37936 - 


-0.30406 


-1.43930 


-2.25393 


-2.71777 


-2.32827 


zp 


10.81108 


4.91586 


l. 75643 - 


-0.10770 


-1.33218 


-2.19131 


-2.67019 


-2.24523 


zv 


7.24295 


3.60651 


l. 97904 


1.21474 


0.82550 


0.63109 


0.69676 


1.55848 


ZA 


6.48301 


3.23151 


l. 79050 


1.H656 


0.77195 


0.59978 


0.67297 


1.51696 



Table 16: S and zx versus smearing parameter for 3 APE clover fermions with csw = l- 



We have also performed a non-perturbative test with csw — 1 clover fermions on our coarsest 
lattice, (3 = 5.846. We find that — am crit decreases monotonically in the range 0<«ape<0.6. 

With ni te r — » oo one expects in perturbation theory that S and [zp — zs)/2 = zy — za tend 
to zero. We checked this explicitly, with details given in Tab.lTTl The approach seems to be 
monotonic in n iter ; we do not observe any oscillations. 



a A ™ = 0.6 


0APE 


1 APE 


10 APE 


100 APE 


1000 APE 


10000 APE 


3 


31.98644 


4.90876 


0.06523 


0.00063 


0.00001 


< 0.0000001 


zs 


19.30995 


4.11106 


-5.94036 


-13.18247 


-20.12297 


-27.03399 


zp 


22.38259 


4.80364 


-5.93562 


-13.18246 


-20.12297 


-27.03399 


zv 


15.32907 


3.31243 


0.16719 


0.01296 


0.00125 


0.00013 


za 


13.79274 


2.96614 


0.16482 


0.01296 


0.00125 


0.00013 



Table 17: 5* and zx versus iteration number for a =0.6 clover fermions with csw = l- 
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